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Abstract 



^ I The quantization scheme is suggested for a spatially inhomogeneous 1+1 

On ■ Bianchi I model. The scheme consists in quantization of the equations of 

^ . motion and gives the operator (so-called quasi-Heisenberg) equations describ- 

^ I ing an explicit evolution of a system. Some particular gauge suitable for 

quantization is proposed. The Wheeler-DeWitt equation is considered in the 
O ■ vicinity of zero scale factor and it is used to construct a space, where the 

^ I quasi-Heisenberg operators act. Spatial discretization as a UV regularization 

procedure is suggested for the equations of motion. 

X. 

b 1 Introduction 



Spatially homogeneous minisuperspace models [TH3] are often used as a testbed for 
the quantum gravity |4]-|9]. Inhomogeneous 1+1 Bianchi I model belongs to the 
so-called midisuperspace models [10] and has more rich properties, in particular, it 
admits an existence of gravitational waves. This model can be reduced to the Gowdy 
one [11] for which the solution of the Wheeler-DeWitt equation has been obtained in 
a closed form [T211I3]. However, the solution of the Wheeler-DeWitt equation does 
not resolve the problem of the gravity quantization completely. An interpretation 
of the Wheeler-DeWitt equation encounters the absence of a variable, which would 
play the role of time, and all approaches to the quantum gravity face, as a rule, this 
challenge. Some approaches regarding the relation of the Wheeler-DeWitt equation 
to dynamics have been suggested. For instance, the approach of Ref. [H] explores 
the notion of "arrival time" from the non-relativistic quantum mechanics to build the 
incoherent histories. Also, there exists a more straightforward approach consisting 
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in quantization of the equations of motion [T514T7] . The result of quantization is the 
so-called quasi-Heisenberg operators. 

Below we apply this approach to the quantization of the Bianchi I model. The 
reason why we investigate the Bianchi I model instead the Gowdy one is that the 
former has a Hamiltonian, which is diagonal on the momentums. Besides, it di- 
vides naturally a spatial geometry into the scale factor and the remaining conformal 
geometry [18], while the Gowdy model suggests other separation. Also, the Wheeler- 
DeWitt equation will be used below since the quasi-Heisenberg operators act in a 
space of solutions of this equation. 



2 Nonuniform Bianchi model 

Let us consider the metric given by the interval 

ds' = e^" (dr]^ - e-^^'dx' - e'^'+'^'^dy' - e^^'-'^'^dz^^ , (1) 

where the functions a,B, V depend on the spatial coordinate x and the conformal 
time rj. The spatially homogeneous metric of such a type has been considered in 
Ref. [19]. 

Substitution of this metric into the Einstein equations allows obtaining a set of 
five independent equations. Two of them are the Hamiltonian and the momentum 
constraints 



2 \6 ' 
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+ ld,,a + hd^Bf + ]-d^,^B + \d,ad.,B + ]-{d.,Vf \ = 0, (2) 

V = e^" (^~d,aa' + d^BB' + '^d^a B' + d^VV + U^B' + U^a'^ = 0, (3) 

and the rest are the equations of motion 

a" - e- {d..a + ^d.af + ^d^Vf + l^d^Bf + \d..B + Ad.Bd.a 

+ V'^ + = 0, (4) 
B" + ^e^^ [d^^B + 2{d,Bf + Q{d^Vf - 2{d,af + 2d,ad^B + 2d,,a) 

+2B'a' = 0, (5) 
V" + 2V'a' - {d,,V + 2d,Vd,a + 4d,Bd,V) = 0. (6) 

The full Hamiltonian H = J l-Ldx has to be zero during an evolution of system. 

The relevant Hamiltonian and the momentum constraints written in the terms of 
momentums 7ry(x) = = e^^V"', tib{.x) = j§f^ = e^"B' and = = 

e^"a' take the form 



2 



V = -padxa + TTsdxB + TivdxV + "^dxTiB + \dxPa- 
Using the Poisson brackets 



^^(x') 5F{x) 5G{x') 6F{x) 6G{x') 



one can obtain the constraint algebra: 

{n{x), nix')} = (P(a;)e^^(") + P(x')e^^("'))5'(x' - x), 
{V{x),V{x)} = {V{x) + V{x))6'{x' - x), 

{'Hix),V{x')} = '^{nix) + 'H{x'))6'{x' - x) - ^'H\x)6{x' - x). (10) 

It is also possible to find the evolution of constraints by calculation of their 
Poisson brackets with the Hamiltonian H: 



d,V{v,x) = {H,Viv,x)} = ^dxnv,x), (11) 



d.n = dx (e^^P) . (12) 



3 Quantization 



The quantization procedure consists in the formulation of initial conditions for the 
quasi-Heisenberg operators. Thereafter it is permissible for operators to evolve in 
accordance with the equations of motion, considered as the operator equations. 
To determine the initial commutation relations, we will use the Dirac quantization 
procedure [201422] . Besides the constraints, an additional gauge condition is needed. 
Let us take the following gauge at the initial moment of time 

A = a — ao, (13) 

B = dxTTB, (14) 

where ao should be tended to — oo. 

For quantization by means of the Dirac brackets [20], one has to calculate the 
matrix Mij{x,x') = {^i{x), $j(x')}, where a set of constraints is 

<I>, = {'H,V,A,B). 
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At the shell of constraints ^i{x) = 0, the matrix Mij{x, x') in the vicinity ao 
— oo has the form 



M{x,x') 



( 



\ 







Pol (x)8{x — x') 

exp(2ao) 








-T^B^"{X - X') 







Pa {x)ll{x-x') 

exp(2ao) 
-^^^ T^B^"{X-X^) 













/ 



• (15) 



Due to antisymmetry of the Poison brackets, Mij{x, x') obeys the identity Mij{x, x') 
—Mji(x', x). The inverse matrix satisfying J Mij(x, x")Mj^{x" , x')dx" = 6ikS{x—x') 
is given by 



M-\x,x') 



( 



\ 






Paix) 









Pa (a;) 3pa(a;)7rs 
g l^{x-x') 



^2an S{x,'-x) !^{x' -x) 



3pa{x')lTB 













where 6{x) is an antiderivative of the Dirac delta-function: 6'{x) = S{x), and A{x) 
is an antiderivative of 6{x): A'{x) = 9{x). It is assumed, that the antiderivatives 
have the following symmetry properties: 0{—x) = —9{x), A{—x) = A{x). 

For the Dirac quantization, one needs to calculate the Dirac brackets which give 
the commutation relations for the corresponding operators at an initial moment of 
time after multiplication by —i. 

Calculation of the Dirac brackets 

{G(x),F(x')}n^{G{x),F(x')} 

~H I {G{x),^i{x'']M;^{x'\x'''){^j{x'''),F{x')]dx''dx''' 



leads to 



{B{x),V{x'))d 



{t:v{x),V{x')]d = Kx-x'), 
{t^b{x),B{x')}d = Q, 

{'Ka{x),a{x')}D = 0, 

{nB{x),Vix')}D = 0, 

{7rBix),Pa{x')}D = 0, 
■Kvix) 



Spaix) TTb 



5{x — x'), 



{B{x),Pa{x')}D = {d-^'Pa{x')e{x' - x) + Pa{x')6{x' - x)) , 



{nB{x),a{x')}D = 0, 
{n^(x),V{x')}n = ^^S{x'-x). 



(16) 



From the foregoing it follows that a and ttb are initially c- numbers in the gauge 
considered. 
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In fact we use some time-dependent gauge, which is known only at an initial 
moment of time, and it is permissable for the commutation relations in the model 
under consideration to evolve. 

Operator realization of the commutation relations at the initial moment of time 
(corresponding to — )■ —oo) may be written as 

... . S 



5V{x) 



B{x) = / S{Tfvix')d^>V{x'))dx' + -pa{,x) 

where the symbol S denotes symmetrization of the noncommutative operators, i.e. 
S{AB) = IIAb + bA) or S{AbC) = HAbC + bAC + ACB + ...). 

Thus the equation of motion (jl]) , ([S]) , dHD should be considered as the operator 
equations of motion with the initial conditions 

y(0,x) = V{x), V'{0,x) = -16-'^"° 



6V{x)' 



3 V SV^{x) 



S'(0,x) =e-2"«7rB, «(0,x) = «o, ^) = "^^7^ + (17) 

where vr^ is some constant and should be tended to — oo. 

Let us suggest a space, where these operators act. Let's consider the Wheeler- 
DeWitt equation in the vicinity of a — — oo 

+ 7r|U[a,y]=o, (18) 



6a{x) 6^V{x) 

where we take into account that vr^ is some constant. 
Solution of Eq. ( [T8|) is of the form of the wave packet 

^[a,V] = J C[iTv] e/(— {-V-l+-?'(-)+-v(x)v^{.))<ix ^^^^ ^^g^ 

Vtiv denotes functional integration over 7iv{x). 

A mean value of an arbitrary operator can be evaluated as [T5HT7] 

< ^\A[a, -i^T, V]\^>=i [( ^*[a, V]D'/^AD-'/^-^^a, V] 
oV J \ oa{x) 
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'^*[a,V]] ID-'/^AlD^/^^[a,V]]vV 



6a{x) ' / ' / 



(20) 



a{x)=ao—^—oo 
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where D{x) = — ^y^^-^ + ttb and W denotes functional integration over V{x). 

In many cases, it is more convenient to use the momentum representation 7rv{x) = 
7ry(x), V{x) = i ^^^(3..) , where the wave function ijj is 

ip[a,7iv] = C['}Tv]exp (-i Oi{x) \J 'Ky{x) + 7r|da; j . (21) 



Then, a mean value of an operator becomes 
< i)\A[a,'nv{,x),i- — — |^ >= 



(22) 



a(a;)=Q!o~^~oo 



Thus, one has an exact quantization scheme consisting of the Wheeler-DeWitt 
equation in the vicinity of small scale factor f[T^ . the operator initial conditions 
(fT7|) for the equations of motion and the expressions fl2U]) . fl22]) for calculation of the 
mean values of operators. 

4 Discretization of the operator equations 

At least two generations of physicists can not overcome divergencies arising in rig- 
orous operator formulation of the ordinary QFT in the 4-dimensional Minkowsky 
space, although some success for the 1+1 and 1+3 dimensional models has been 
reached [23]. Here we consider a discretization as a method of regularization of the 
functional operator equations, which eliminates the infinite quantities. An alter- 
native is to use the discretized action initially, however, we do not consider such a 
possibility here. The discretization consists in choosing of some spatial box of length 
L and granulation of it by points Xi separated by distance Aa;. Periodicity condition 
is implied aX, xi = xjs;. 

Continuous oscillators of \^-field should be replaced by the discrete momentums 
7ry(a;) — )■ T^vj/^ Ax, V{x) — )■ Vj/ V Ax [21], where vryj and Vj posses an ordinary 
commutation relation [vry^, Kn] = —i^nm including the Kronecker symbol 5mn- Re- 
ally in this case for 7ry(x„) and V{xm), we have [7ry(x„), V"(xm)] = —i5nml 
that turns into the Dirac delta-function bnmj Ax — > 5(x„ — Xm) in the limit Ax — )■ 0. 
However, it is more convenient to use straightforwardly the quantities vryj = 7ry(xj), 
Vj = V{xj) for which the commutation relations [7ry„, V^] = —i6nm/Ax are satis- 
fied. These commutation relations can be realized by the operators = 
7ry„ = -i^^ or vry^ = 7ry„, K = ^g^. 

The discrete Wheeler-DeWitt equation in the vicinity of aj = ao ^ —00 has the 
following form in the momentum representation: 

+ Tiy^ + 7r| 4){ai, . . . aN, t^vi, ■ ■ ■ '"'vn) = 0, (23) 



(Ax)2 d^a. 
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with a solution 



ilj{ai . . . On, TTyi, . . . , hvn) 



N 



C{7vvu , TTvTv) exp ^-iAx ^ c^jyT^vj + T^Bj ) (24) 

where Ax is the discretization length. 

A mean value of an arbitrary operator can be evaluated as 

i d i 

< V'|A(ai...aAr,7rv'i,...7ryAf, — > = 

Ax OTTvi AxonvN 

C{tTvi ■ ■ ■'JTvN)dTTvi ■ . . dlTvN -(25) 

ai...ajv=c«0— 

For the equations of motion, one has to take some operator ordering and then to 
come to its discrete version. A question arises about the conservation of constraints 
during evolution. Constraints can be violated by both noncommutativity of opera- 
tors and discretization of the equations of motion. The last question is analogous 
to that regarding the energy conservation in a system of the decretized magnetic 
hydrodynamic equations [25] . Special discretization schemes were suggested for the 
magnetic hydrodynamic equations, which conserve the energy |25j . 

Let us propose a solution of both problems. Let's consider the discretized oper- 
ator Hamiltonian, the momentum constraints and, besides, choose the symmetrical 
ordering of operators: 



^ \ \ 3Ax ^ Ax ^ 3 Ax ^ 

+ Ax 3Ax + 3Ai ))' ^ ' 

% = S (ie-- (-af + Bf + V-) + e^'*'"' [^i^) ' 



3(Ax)2 6 \ Ax / 3(Ax) 



2n 



+^ "'""til^"""'' ^H^' '''' 

Here, we need more careful definition of the symmetrization [26l|27]. Let's there 
is an arbitrary function /(xi, X2, ■ ■ ■ x„) of n— variables. Then, one can define a 
formal Fourier transform 

/(Ci, C2 . . . Cn) = 7^ / /(xi, X2, . . . x„)e-^(^^^^+^^^-^"^")rfxi . . . dx^. (28) 



7 



A symmetrized function of noncommutative operators Ai, . . . A„ is defined as 

S{f{A„ A,... i„)) = I /(Ci, C2, . . . Cn)e"^^^^'+^'<'-^-<-^dC, . . . dCn. (29) 



Our idea is to use the discretized version of Eqs. (jTTj) . ( 1121) as the equations of 
motion: 



d,H, = ' ' > (30) 

Using the formula for differentiation of a symmetrized function [27]: 

pifiMt), Mt) ■ ■ ■ Mt))) = s(^^ d,f{A,{t), A,{t) . . . i„(t)) j , (32) 

(here djf{xi, X2 ■ ■ ■ Xn) denotes a partial derivative of a function / over the j— argument) 
allows calculating the time derivatives in the left hand side of Eqs. ( l30|3Ti) and 
rewriting Eqs. ( 1301I3T1) in the form of 

s (^e^^^ (^^b; + v;v;' - + ^e^^^ (jb] + tb^, + 2b,_, + 3v- 

-24(2 + 74+1 - 4dj + 4dj+i) + 24+1(1 - 4dj + A&j+i) - 6VjVj+i 
+31^+1 — 4aj + ctj + 2dj+i — 2ajaj+i + ctj^i + 2dj_i 1 I 2Bj + 1 

+«; + - «?) + ^^^'"^^ ((-2 + 74 - 74+1 + 4A, - 4a,+04 

+(1 - 74 + 74+1 - 4a,- + 4d,-+i)4+i + 4-1 + 3v;-v:' - 3v:-+iV:' 



(1 - 74 + 74+1 - 4a,- + 4d,-+i)4-+i + 4-1 + 3^i^, 
3VjVj_^_i + 3Vj-+iVj+]^ — 2d^- + ABja'j — iBj^ia'j + ajaj — aj+id^ + d^+j^ 

— 44dj+i + 44+ id^+i — djd^-4.^ + dj+id^_^i + d^_i 
_ + 4+16^^^+^ - - 4e'^- 



2Ai • '^^) 



'J 



+3v;iv;^, + 2(-v;- + r,+i)d;.) + d;d;+i + 4(34+1 + (-4 - 64 + 64+1 - 4«,- 

+4d,-+i)d;. + 2d;.+i) + (-1 - 34 + 3fi,+i)4' - 2d,-4' + 2d,-+i4' + 4>i 

+3{-v, + v,+,)v;' - a'; + (d, - d,-+i)d;' + a';^,^ j = 4+1 - 4, (34) 
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where l-ij and Vj are given by fl26l) . fl271) . Moreover, one can consider Eqs. fl33l) . 
(IMj) without the right hand side as the equations of motion, since i-Lj, Vj equal to 
zero initially at ao ~ ~oo in accordance with the initial conditions. One must be 
sure only that the time derivatives V'^ equal to zero, as well. 

As the third equation, the discretized version of Eq. (jS]) can be taken 



[Ax)'^ (Ax)"' 




Thus we have three operator equations ( !33|) . ( !34j) . ( 1351) (or equivalently ( !MI) .(l35l) 
without the right-hand side) which for commuting quantities and in the continuous 
limit Ax — )■ are fully equivalent to the classical equations ([5]), (j6]). 

The equations should be solved with the following initial condition 



4(0) 




d 



A;. = e-^-Jvr^^.+7r|. (36) 



In the vicinity of small scale factors, the time derivatives are much larger than 
the spatial ones so that the operator equations take a simple form 



+ + V;'2 + 5f = 0, (37) 

B] + 1B)6lj = 0, (38) 

V;' + 2V^a'^ = 0, (39) 

with the constraint —af + Bf + V-"^ = 0. All quantities, i.e. a'^B'^V- in these 
equations commute with each other. 

The solution of Eqs. (!37|) . (137|) . (138!) with the initial conditions (136|) is given as 

V,{V) = Vm + Z"^^' In (l + 2e-^-' J^^lT^r/ 



BM = 4(0) + ^ In (l + 2e-2"y vrl + TT^.^. r^) , 



= «o + ^ In (l + 2e-2"°^/^4T4~r7 



Although the evolution in the vicinity of a ~ ao ~ "Oo is relatively simple, the 
evolution governed by the general equations (!33|) . (!34|) . ( 135|) . when the fields begin 
to oscillate is very complicated and needs the numerical investigation that is beyond 
the aims of the present article. 
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5 Discussion and Conclusion 

It should be said about a choice of a state of system considered. We do not iden- 
tify the state as an "initial" because there is only the state C^tcvi, ,''^vn), 

which describes all evolution of system and allows calculating the mean values of 
the quasi-Heisenberg operators. This state is not related to the notion of "vacuum" 
since there is no field oscillators in the limit a — > — cxd. In Ref. pjj solution of 
the Wheeler-DeWitt equation for the Gowdy model is investigated. Although the 
another quantities were introduced in [TD] instead of logarithm of scale factor a and 
-B-field, namely 



the asymptotic of the Wheeler-DeWitt equation in the vicinity T — oo has the 
form containing only the momentums by analogy with Eq. (IT5]) and admits solutions 
of the plane wave type. Asymptotic of the Wheeler-DeWitt equation in the vicinity 
T — >^ oo is of the oscillator type |T9]. Quasiclassical treatment of the Wheeler- 
DeWitt equation with regard to an evolution of universe allows an interpretation as 
a scattering problem [19] i.e. transition of the packet of plane waves at T — )■ — oo 
to the number of gravitons (i.e exited oscillators) at T — )• oo. It should be noted 
that there were no the state at T — !■ — oo which give no gravitons at T — )■ oo. The 
later work [T3| about a quantization of the Gowdy model was in some sense a step 
back because it considered graviton creation from vacuum in a style of Refs. [281431] . 
It was suggested that at some time Tq where the field oscillators already have 
formed there is no gravitons and later the gravitons appear from vacuum during 
an evolution. However, our opinion is that there is no physical background for the 
existence of such a time Tq where the universe was empty. 

In our evolutionary picture, we have also the state describing evolution of uni- 
verse as a packet of the "plane waves" . In a future this will give some gravitons under 
the vacuum. More exactly if one calculated the correlators < V{r],x)V{ri,x') > he 
will find that they are analogous to the correlators of the QFT corresponding to 
some gravitons under vacuum. There is no wave packet which gives a pure vacuum 
of universe in the future, i.e. appearance of matter is inevitable in this model but 
the matter is not created from vacuum, because at any time the universe is not 
empty. 

As an example of the initial state one may take 



where Cat is the normalization constant. This state implies that the momentums 
are random and independent in spatial points. This state does not look similar to a 
vacuum state of the QFT, as the fields (and their momentums) in a QFT vacuum 
state are highly correlated in the nearest spatial points. 

As an example of the mean value calculation, one can take the evolution of a 



T = B + a, 
A = 6{B-a) 



(40) 




(41) 
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"Hubble constant" of system at small rj: 



1 ddj 1 ddj , „ .doLj h^l'^ /, o/. 7 2n /.on 

< af Af ^^P<-"') Af >= v7it'(V4.3/4>.J). (42) 

where f/(a, 6, is the confluent hypergeometric function. 

For this state, the universe is expanded uniformly in a mean. It would be inter- 
esting to calculate the evolution of correlators < yj{j])Vn{j]) >■ Initially the fleld 
Vj is uncorrelated for different j: < Vj{ri)Vn{ri) >~ 6jn but then some correlation 
should arise and some analogy with the vacuum of QFT in this sense should appear 
so that the flelds become highly correlated in the nearest spatial points. 

The question arises: why do the commutation relations related to the fleld V 
have the usual form initially but do become to be broken after an evolution and, 
probably, look differently (the difference have to be related to the Hubble constant) 
at the present time? The answer is that it is a consequence of the chosen initial 
gauge. The gauge can be different, but the gauge considered is simplest and allows 
obtaining the commutation relations explicitly at the initial moment of time. As a 
matter of fact, the gauge choice remains to be complex issue. 

To summarize, failure of the QFT in flat spacetime to deal with such an inher- 
ently non-linear theory as gravity and existence of "problem of time" insist on an 
invention of some new quantization procedures. The considered quasi-Heisenberg 
quantization scheme may provide a calculational framework for the investigation of 
the quantum evolution of a system. The goal of further investigation may be the 
vacuum energy problem, more exactly, its possible zero value in the quantization 
scheme considered, that may result from compensation of the zero point fluctuations 
of gravitational waves by the quantum fluctuations of the scale factor. Thereby, the 
fluctuations do not contribute to the mean evolution of a system. It should be 
noted that this will be purely quantum effect, since it is absent in classics [32j. An 
additional issue is calculation of the fleld correlators in order to determine their 
correspondence to the correlators of ordinary QFT in late times. 
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